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ABSTRACT 

We analyze the evolution of 42 spiral galaxies in the Spitzer Infrared Nearby Galaxies Survey 
(SINGS). We make use of ultraviolet (UV), optical and near- infrared radial profiles, corrected for 
internal extinction using the total-infrared to UV ratio, to probe the emission of stellar populations 
of different ages as a function of galactocentric distance. We fit these radial profiles with models 
that describe the chemical and spectro-photometric evolution of spiral disks within a self-consistent 
framework. These backward evolutionary models succesfuUy reproduce the multi-wavelength profiles 
of our galaxies, except the UV profiles of some early-type disks for which the models seem to retain 
too much gas. From the model fitting we infer the maximum circular velocity of the rotation curve 
Vc and the dimensionless spin parameter A. The values of Vc are in good agreement with the 
velocities measured in HI rotation curves. Even though our sample is not volume-limited, the resulting 
distribution of A is close to the lognormal function obtained in cosmological A^-body simulations, 
peaking at A ^ 0.03 regardless of the total halo mass. We do not find any evident trend between 
A and Hubble type, besides an increase in the scatter for the latest types. According to the model, 
galaxies evolve along a roughly constant mass-size relation, increasing their scale-lengths as they 
become more massive. The radial scale-length of most disks in our sample seems to have increased 
at a rate of 0.05-0.06 kpc Gyr~^, although the same cannot be said of a volume-limited sample. In 
relative terms, the scale-length has grown by 20-25% since z = 1 and, unlike the former figure, we 
argue that this relative growth rate can be indeed representative of a complete galaxy sample. 
Subject headings: galaxies: abundances — galaxies: evolution — galaxies: photometry — galaxies: 
spiral 



1. INTRODUCTION 

Unveiling the details that govern the formation of disk 
galaxies is paramount for our understanding of the evo- 
lution of the universe as a whole. In the currently ac- 
cepted paradigm of galaxy formation, rotating proto- 
galactic clouds collapse within the gravitational wells of 
dark matter haloes. Gas cools via radiative processes 
I and, if it keeps enough angular momentum, a rotationally 
supported gaseous disk will eventually form (Fall & Efs- 

■ tathiou 1980; White & Frenk 1991; Mo et al. 1998). Dark 
matter haloes themselves grow from primordial density 
fluctuations, and they are supposed to merge and evolve 

' according to the A Cold Dark Matter model (ACDM; 

■ Springel et al. 2005; Spergel et al. 2007). 
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Because gas takes longer to settle onto the disk in the 
outer parts, given its larger angular momentum and the 
longer gravitational collapse time, star formation should 
proceed on longer timescales in the outskirts of disks than 
in the inner regions. Therefore, a natural consequence of 
such a scenario is that disk galaxies should be assembled 
from inside out (Samland & Gerhard 2003). In partic- 
ular, the radial scale-length of exponential disks is ex- 
pected to increase with time (Brook et al. 2006; Brooks 
et al. 2010). In principle, the mass and size evolution 
of galaxies can be probed with observations at different 
redshifts (see, e.g. Trujillo et al. 2004, 2006; Bardcn et 
al. 2005 and references therein), provided that one can 
properly deal with cosmological and selection effects. 

In a way complementary to this look-back approach, 
galactic archaeology in the local universe has also proven 
useful to infer the past evolution of galaxies. Each partic- 
ular scenario of galactic evolution should have left charac- 
teristic imprints in the radial variation of the properties 
of stars, gas and dust in present-day galaxies. If galaxies 
do indeed grow from inside out, stars should be younger 
on average in the outer parts, leading to radial color gra- 
dients such as those we actually observe (de Jong 1996; 
Beh & de Jong 2000; MacArthur et al. 2004; Taylor et 
al. 2005; Mufioz-Mateos et al. 2007). Measuring age gra- 
dients in disks from color profiles is not straightforward, 
since the radial decrease in the internal extinction and 
metallicity also conspire to yield bluer colors at larger 
radius. 

More recently, color-magnitude diagrams of resolved 
stellar fields in nearby galaxies have also favored an 
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inside-out scenario of galactic growth (Gogarten et al. 
2010; Barker et al. 2010). There is also ample observa- 
tional evidence that chemical abundances decrease with 
galactocentric distance, both in the Milky Way (MW; 
Shaver et al. 1983; Smartt & RoUeston 1997) and in 
external galaxies (Zaritsky et al. 1994; van Zee et al. 
1998; Pilyugin et al. 2004; Moustakas & Kennicutt 2006; 
Moustakas et al. 2010). Multi-zone chemical evolution 
models based on the inside-out scenario usually invoke 
a radially-increasing timescale of gas infall to reproduce 
these metallicity gradients (Matteucci & Francois 1989; 
MoUa, Ferrini & Diaz 1996; Prantzos & Boissier 2000; 
Chiappini et al. 2001; Carigi et al. 2005). 

In this paper, the third in a series devoted to the spa- 
tial distribution of stars, gas and dust in nearby galaxies, 
we will test the predictions of the multi-zone model of 
Boissier & Prantzos (1999, 2000; BP99 and BPOO here- 
after). This model describes the chemical and spectro- 
photometric evolution of spiral disks in a self-consistent 
framework, taking into account the radially-varying gas 
infall rate, a physically motivated star formation law and 
a full treatment of chemical evolution. The model is first 
calibrated to reproduce observables of the MW (BP99) 
and then extended to other disk-like galaxies through 
scaling laws resulting from the ACDM model (BPOO). 
Using only two free parameters, the maximum rotational 
velocity of the rotation curve Vc and the dimensionless 
spin parameter A, the model is able to predict radial 
profiles of several quantities, including multi-wavelength 
photometry, metallicity and gas density among others. 

In Mufioz-Mateos et al. (2009a, Paper I hereafter) we 
derived multi- wavelength profiles from the far-ultraviolet 
(FUV) to the far-infrared (FIR) for the galaxies in the 
Spitzer Infrared Nearby Galaxies Survey (SINGS; Ken- 
nicutt et al. 2003), which comprises 75 objects represen- 
tative of the typical galaxy population in the Local Uni- 
verse. In Muiioz-Mateos et al. (2009b, Paper II hereafter) 
we centered our attention on the radial variation of sev- 
eral physical properties of dust, such as the internal ex- 
tinction, the dust mass surface density, the abundance of 
polycyclic aromatic hydrocarbons (PAHs) and the dust- 
to-gas ratio. In the present paper we have focused on a 
subsample of 42 spiral galaxies within the SINGS sample. 
We have combined the UV, optical and near-IR profiles 
measured in Paper I with the extinction profiles obtained 
in Paper II, in order to recover the intrinsic emission of 
stars of different ages across the galactic disks. These 
stellar, extinction-free profiles have been fitted with the 
models of BPOO, thus testing the ability of the models to 
simultaneously reproduce the multi-band profiles, while 
at the same time inferring the circular velocities and spin 
parameters of each object. 

The purpose of the present paper is thus twofold: first, 
we will verify whether the models are able to reproduce 
the present-day profiles of nearby disks. Second, we will 
indirectly obtain the values of Vc and A for each galaxy. 
The spin parameter is particularly important in cosmo- 
logical studies. While Vc can be easily determined from 
rotation curves, A is not a directly measurable quantity. 
Previous studies (Syer et al. 1999; Hernandez et al. 2007; 
Cervantes-Sodi et al. 2008) have shown that A can be em- 
pirically estimated from a combination of observed galac- 
tic properties such as the disk scale-length, provided that 
some ACDM-based assumptions are made. When apply- 



ing this methodology to large optical datasets of nearby 
galaxies, these authors found an excellent agreement be- 
tween the empiricial distribution of A values and the one 
obtained in iV-body simulations of hierarchical cluster- 
ing. However, optical measurements of disk scale-lengths 
might be biased by radial variations of the mass-to-light 
ratio or the internal extinction. The BPOO models incor- 
porate the radial variation of gas infall (inside-out forma- 
tion) and of the star formation law, thereby accounting 
for wavelength variations in the disk scale-length in a nat- 
ural way. Moreover, the surface brightness profiles we use 
to constrain the models are corrected for internal extinc- 
tion using robust methods (see Paper II). Besides, apart 
from UV and optical profiles, we incorporate in our anal- 
ysis near-IR ones, which are less sensitive to mass-to- light 
variations and dust attenuation. This extra wavelength 
coverage obviously comes at the expense of using a much 
more reduced sample than in the aforementioned studies. 
Therefore, our analysis cannot reach the same levels of 
statistical completeness, but it should serve nevertheless 
as a robust foundation for future works. 

Features like bulges, bars and radial mass flows are not 
considered in the BPOO models. Accounting for all these 
effects require A'^-body simulations, which given their 
complexity are usually limited to a handful of objects. 
Therefore, despite their somewhat simplified underlying 
assumptions, models such as the BPOO ones allow simu- 
lating large grids of galaxies that cover a wide range of 
properties, and are therefore better suited for our pur- 
poses. 

This paper is organized as follows. In Section [5] we pro- 
vide a brief summary of the inner working of the BPOO 
models. In Section [3] we describe the main properties 
of our subsample of SINGS spirals, and summarize the 
procedure used in Paper I to derive the surface bright- 
ness profiles. In Section |4] we explain how the multi- 
wavelength profiles were corrected for the effects of in- 
ternal extinction. Section [S] deals with the details of 
the fitting procedure, and in Section |6] we present the 
main results of our analysis. The main conclusions of 
this work are summarized in Section [T] Finally, the two- 
dimensional distribution of values resulting from the 
model fitting are compiled in Appendix |^ 

2. DESCRIPTION OF THE MODELS 

In this section we broadly outline the main in- 
gredients and underlying assumptions of the chemo- 
spectrophotometric models used to fit the multi- 
wavelength profiles of the SINGS galaxies. The reader 
is referred to BP99 and BPOO for a more in-depth de- 
scription of the physical details of the models. Briefly, 
an initial model was first developed and calibrated to re- 
produce several observed properties of the MW (BP99). 
This model was then generalized to other spiral disks 
of different sizes and masses by means of several scaling 
laws deduced from the ACDM scenario of disk formation 
(BPOO). 

2.1. The Milky Way model 

The Milky Way disk is modeled as several concentric 
rings, which are progressively built up by accretion of 
primordial gas from the halo. In the model, these an- 
nuli evolve independently from one another, in the sense 
that no radial mass flows are allowed. Such flows can 
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actually take place in real galaxies as a result, for in- 
stance, of the presence of bars (Sellwood & Wilkinson 
1993), redistribution of angular momentum due to viscos- 
ity (Yoshii & Sommer-Larsen 1989; Ferguson & Clarke 
2001) and radial stellar migration (Roskar et al. 2008; 
Martinez-Serrano et al. 2009; Sanchez-Blazquez et al. 
2009). In particular, stellar migration has been proposed 
as a likely mechanism to explain the observed U-shaped 
color profiles in galaxies (Azzollini et al. 2008a; Bakos 
et al. 2008). Moreover, gas outflows are not included in 
the models, while observations suggest that they play a 
role in the chemical evolution of low-mass galaxies (Gar- 
nett 2002). Nevertheless, this is likely not a concern in 
our analysis, since most of our disks are large and mas- 
sive enough. Finally, the model does not include the 
bulge nor it does differentiate between the thin and thick 
disk. Implementing all these phenomena in an analytic 
way is not straightforward, and it would introduce many 
additional free parameters whose values might be diffi- 
cult to constrain. Despite the simplifying assumption of 
independently-evolving rings, the model is still successful 
at reproducing the radial structure of the Milky Way. 

The star formation rate (SFR) surface density at each 
radius r and time i, I]sFR(ij'')j depends on the local gas 
density Sg(i, r) following a Schmidt law modulated by a 
dynamical term: 

SsFR(t,r) =aI]g(t,r)"F(r)r^i (1) 

Here V{r) is the rotational velocity at radius r. The 
term V{r)r~^ is intended to mimic the conversion of gas 
into stars by the periodic passage of spiral density waves 
(Wyse & Silk 1989); it can be also seen as the inverse of a 
dynamical timescale. Originally, a was fixed in order to 
reproduce the local gas fraction in the solar neighbour- 
hood at T = 13.5 Gyr, and n = 1.5 was chosen to re- 
produce radial trends. Later on, in Boissier et al. (2003) 
the star formation law was empirically determined using 
a sample of nearby galaxies. Here we adopt the values 
of a = 0.0026^ and n = 1.48 found in that work. Note 
that these values are very close to the ones originally 
adopted in BP99 and BPOO [a = 0.00364 and n = 1.5). 

The gas infall rate f{r,t) decreases exponentially with 
time: 

/(t,r)=A(r)e-*/^M (2) 

The timescale of gas accretion T(r) is assumed to in- 
crease with radius, from 1 Gyr at r = 1 kpc to 15 Gyr 
at r = 17 kpc. This allows reproducing the inside- 
out formation of disks, since gas settles onto the disk 
on longer timescales in the outer regions due to having 
larger angular momentum. At r = 8 kpc r is set equal 
to 7 Gyr to reproduce the metallicity distribution of G- 
dwarf stars in the solar neighbourhood. The normalizing 
factor A{r) can be deduced by integrating the infall rate 
until T = 13.5 Gyr, and then matching the result to the 
current stellar mass profile of the disk: 

r f{t,r)dt = Esoiar ncigh.e-^''-'*'^P^^/''^^ = Soce"''/^'^'^ 

"'0 



The units of a are such that EspR(r) is measured in 
Mq pc~^ Gy r~^ , Sgas (r) in Mq pc~^, r in kpc and V{r) in km s~^. 
See also Fig. \U\ 



The subscript G refers to the parameters of our galaxy. 
According to observations of the Milky Way, the radial 
scale-length is fixed to i?dG = 2.6 kpc, and the central 
mass density (extrapolated from the one in the solar 
neighborhood) is set to Sqg = 1150 Mgpc'^ (see BP99 
for references). 

The distribution of stars for a given SFR follows a user- 
specified initial mass function (IMF) . Even though we are 
explicitly assuming the existence of a universal IMF, this 
might not be necessarily the case (see Bastian et al. 2010 
for a review on the subject). Therefore, in this work 
we will compare the results obtained with the IMFs of 
Kroupa et al. (1993; K93 hereafter) and Kroupa (2001; 
KOI hereafter). The K93 IMF was used in the original 
models, but here we are also interested in analyzing the 
results yielded by a more recent version of the IMF. The 
optical and near-IR fluxes, as well as the gas quantities, 
change by less than 20% between these two IMFs, which 
justifies not to totally recalibrate the model. It is the 
UV fiuxes and metallicities that vary significantly, given 
the different content in high-mass stars of both IMFs 
(see section [6.2.3p . Of course, there exist several other 
parameterizations of the IMF which are widely used in 
extragalactic studies, but the comparison between the 
K93 and KOI IMFs presented here should suffice for the 
purpose of showing the impact of varying the relative 
amount of high- and low-mass stars. 

Stars of different masses enrich the ISM with varying 
amounts of different elements; in this regard, the model 
does not assume the instantaneous recycling approxi- 
mation (Tinsley 1980), according to which stars more 
massive than 1 Mq die instantly, whereas less massive 
ones live forever. On the contrary, the model takes into 
account the finite lifetimes of stars of different masses 
when computing the chemical evolution within each ring. 
Moreover, the properties of each new generation of stars 
(lifetimes, stellar yields, evolutionary tracks and spectra) 
depend on the local metallicity at the corresponding ra- 
dius and time of formation (see BP99). The spectrum 
of a given ring at time t can be then computed as the 
sum (both in time and mass) of the individual spectra of 
previously formed stars which are still alive at time t. 

With the assumptions outlined above, BP99 showed 
that their Milky Way model is able to reproduce not only 
observables in the solar neighbourhood, but also radially- 
dependent ones, such as profiles of gas surface density, 
gas-phase oxygen abundance, SFR and supernova rates, 
as well as luminosity profiles at different bands. 

2.2. Extension to other disk-like galaxies 

The previous model for the Milky Way was general- 
ized to other disks in BPOO, by making use of the scaling 
laws derived by Mo et al. (1998) within the ACDM sce- 
nario. In this theoretical framework, galaxy formation 
is usually split into two different processes: the growth 
of non-baryonic dark matter haloes and the assembly 
of baryonic structures within them. Gravitational in- 
stabilities amplify the primordial density fluctuations, 
yielding dark matter clumps that merge and interact 
with each other, acquiring angular torques during the 
process. Meanwhile, baryonic gas cools and condenses 
within these haloes, leading to self-gravitating structures 
that are able to form stars, thus eventually giving rise to 
present-day galaxies. 
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The models of BPOO build on the mathematical for- 
malism of Mo et al. (1998), which establishes that un- 
der certain assumptions the scaling properties of disks 
depend only on two parameters: the maximum circular 
velocity of the rotation curve Vq and the dimensionless 
spin parameter A: 



Vc = [10G'i/(z)M]i/3 



(4) 
(5) 



In the equations above, M, J and E are the total mass, 
angular momentum and energy of the halo, G is the grav- 
itational constant and H(z) is the Hubble parameter at 
the redshift z of halo formation. In order to express the 
properties of disks in terms of Vq and A alone, the fol- 
lowing assumptions need to be made: 

1. The masses of disks Afd are just a few percent of 
those of their corresponding haloes. The precise 
value of this ratio is unclear, but it must conform 
to the baryonic fraction of the universe and the 
efficiency of disk formation. Following Mo et al. 
(1998), the BPOO models assume Md = 0.05Af for 
all disks. 

2. The specific angular momentum of the disk and 
halo are equal (i.e., Jd/Md = J/M). While this 
commonly used assumption is not strictly sup- 
ported by numerical simulations, it is apparently 
required to produce disk sizes that match observa- 
tions. 

3. Variations in the formation time of the disks are ig- 
nored. It is now believed that the thin component 
of disks is assembled at ^ ^ 1 (Brook et al. 2006), 
and its evolution dominates the inside-out growth 
of spirals until z = (Chiappini et al. 1997). How- 
ever, disks might contain stellar populations that 
formed much earlier. The concept of 'formation 
time' is thus somehow ill-defined, and the BPOO 
models simply assume that all disks started form- 
ing stars at the same time, having today a fixed 
age of 13.5 Gyr. 

Under these assumptions, BPOO showed that the scale- 
length i?d and central mass density Eq of a given disk can 
be derived from those of the Milky Way by means of their 
relative spins and circular velocities: 



So 



A Vc 
Ag Vcg 

A' ^ 

Ag 



Vc 

VcG 



(6) 
(7) 



For the case of the Milky Way, the BPOO models assume 
that VcG = 220kms"i and Aq = 0.03. Although both 
Vc and A affect the final scale-length of a disk, they do 
it in different ways, as can be seen in Fig. [TJ Larger 
values of Vc yield more extended and massive disks, while 
modifying A alters the scale-length alone. 

Note that when calibrating our model on the MW, we 
are implicitly assuming some sort of homology between 
the evolution of the galaxies we study and the MW. This 
homology is observationally motivated by the fact that 
several properties of galaxies such as metallicity gradients 
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Fig. 1. — Sample surface brightness profiles generated by tfie 
model at 3.6 lira, showing the effect of varying the spin parameter 
(a) and the circular velocity (b). 

(Garnett et al. 1997; Henry & Worthey 1999) or rotation 
curves (Salucci & Persic 1997) become 'universal' once 
normalized by the optical size. Nevertheless, homology 
constitutes a symplifying assumption that may be partly 
responsible of some of the discrepancies that we will show 
later. 

The scaling laws described above affect the way in 
which the SFR and the gas infall time-scale depend on 
galactocentric distance. The final rotation curve of a 
given simulated galaxy is computed as the sum of the 
contributions of the halo and the disk. The resulting 
function V{r) is then used through Eq. [T] to determine 
the radial variation of the SFR. The time-scale for the 
gas infall is parameterized as a function of both the local 
mass surface surface density and the total galaxy's mass, 
in the sense that a deeper gravitational well leads to a 
more rapid infall of gas onto the disk (see Eq. 21 and 
Fig. 3 in BPOO). 

3. THE DATA 

From the original 75 objects of the SINGS sample we 
first exclude all ellipticals, lenticulars and dwarf irregu- 
lars, leaving only those galaxies with morphological types 
1 < T < 9. From the remaining list of galaxies we also 
exclude the following objects: 

1. NGC 2798, a Sa galaxy with a severely distorted 
morphology due to its interaction with the neighbor 
galaxy NGC 2799. 

2. NGC 3190 and NGC 4594 (The Sombrero Galaxy). 
These Sa galaxies are seen almost edge-on, with 
dense dust lanes heavily obscuring part of their 
disks. Besides, their prominent bulges modify the 
ellipticity of the isophotes used to measure their 
surface brightness profiles, which may not be rep- 
resentative of their disk components alone. 
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3. NGC 4631, an edge-on Sd galaxy for which our 
elliptical isophotes probably mix light emitted at 
very different galactocentric distances. 

4. NGC 5474, a Scd galaxy with a disturbed mor- 
phology, probably due to a tidal interaction with 
M 101. In the optical and near-IR, its main disk 
is significantly shifted southwards with respect to 
the bulge. 

Even though NGC 5194 (M51a) is cleary interacting 
with NGC 5195 (M51b), we do not exclude it from our 
analysis, since it still retains it axial symmetry. After ap- 
plying these criteria we are left with 42 disk-like galaxies, 
whose main properties are summarized in Table [1] 

We refer the reader to Paper I for a detailed descrip- 
tion of the imaging dataset employed here. UV images 
were taken with the GALEX space telescope (Martin et 
al. 2005) in the FUV and NUV bands, and belong to the 
sample compiled in the GALEX Atlas of Nearby Galax- 
ies (Gil de Paz et al. 2007). We employed two sets of 
optical images. On one hand, we relied on ugriz images 
from the Sloan Digital Sky Survey DR6 (SDSS; York et 
al. 2000; Adelman-McCarthy et al. 2008) when available. 
For objects without SDSS imaging, we used the origi- 
nal SINGS optical images (Dale et al. 2007), taken with 
the Kitt Peak National Observatory (KPNO) 2.1 ni tele- 
scope and the Cerro Tololo Inter- American Observatory 
(CTIO) 1.5 m telescope. Some of these BVRI optical 
images were affected by zero-point offsets that were cor- 
rected in Paper I. Note that this recalibration procedure 
adds an extra component to the photometric uncertainty 
of the BVRI fluxes that is not present in the SDSS ones. 

Near-IR images in the J, H and Ks bands were com- 
piled from the 2MASS Large Galaxy Atlas (LGA; Jarrett 
et al. 2003). Finally, we also used images at 3.6 /im and 
4.5 /zm, taken with the Infrared Array Camera (IRAC; 
Fazio et al. 2004) onboard Spitzer (Werner et al. 2004). 
Since we are interested in tracing the stellar emission, we 
did not consider the other two IRAC bands at 5.8 /xm and 
8.0 /im, which contain significant emission from hot dust 
and PAHs. However, we did use those bands, together 
with the FIR bands at 24, 70 and 160 fiin from the Multi- 
band Imaging Photometer (MIPS, Rieke et al. 2004) to 
compute the radial variation of internal extinction from 
the TIR/FUV and TIR/NUV ratios (see Paper II and 
Section [31). 

Technical details on how the radial profiles were ob- 
tained are also given in Paper I. Briefly, we used the 
IRAlfl task ELLIPSE to measure the mean surface bright- 
ness along concentric elliptical isophotes, using the same 
sets of ellipses at all bands for each galaxy. The ellipticity 
and position angle were kept fixed and equal to those of 
the /is = 25 mag arcsec"^ isophote from the RC3 catalog 
(de Vaucouleurs et al. 1991), or from the NASA Extra- 
galactic Database (NED) when these parameters were 
not included in the RC3. They are quoted in Table [T] to- 
gether with the central coordinates of the ellipses, which 
were also maintained fixed. We used radial increments of 
6" along the semimajor axis (similar to the FWHM of the 

^ IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



GALEX PSF) up to a final radius at least 1.5 times the 
diameter at /is = 25 mag arcsec"^ (D25). This upper 
limit was increased when significant emission was seen 
beyond that radius (especially in the UV bands). The 
uncertainty of the mean surface brightness at each ra- 
dius comprises the Poisson noise in the source flux and 
the error in the sky level, the latter including both local 
and large-scale background variations (see Paper I) . 

4. INTERNAL EXTINCTION CORRECTION 

Prior to fltting the multi-wavelength profiles, we must 
first correct them for the radial variation of the inter- 
nal attenuatiorPI. In Paper II we computed internal 
attenuation profiles in the FUV and NUV bands inde- 
pendently from the TIR/FUV and TIR/NUV ratios, re- 
spectively. We followed the prescriptions of Cortese et al. 
(2008), which take into account the varying extra dust 
heating due to evolved stellar populations. After ^fuv 
and Anuv have been obtained, the attenuation at other 
wavelengths can be derived after assuming a given ex- 
tinction law and a geometry for the distribution of stars 
and dust. Here we follow the prescriptions of Boselli et 
al. (2003) and adopt a sandwich model, where a thin 
layer of dust is embedded in a thicker layer of stars: 



A,(A) = -2.51og 



l-C(A) 



C(A) 



r(A) sec(i) 



I g— r(A)scc(i) 



-r(A) scc(i) 



(8) 



Here t(A) is the face-on optical depth and i is the incli- 
nation angle. Note that we need not to care about them 
separately as t(A) sec(i) is a joint quantity. The variable 
C(A) denotes the ratio between the thickness of the dust 
and stars layers. Young stars, which dominate the emis- 
sion in the UV range, are likely immersed in a thin dust 
layer. More evolved ones, which emit most of their light 
predominantly in the optical and near-IR bands, migrate 
with time out of the galactic plane, and are thus assumed 
to lie within a thicker layer, partly above and below the 
thin dust layer. Therefore, Boselli et al. (2003) parame- 
terize the dust-to-stars scale-height ratio as a decreasing 
function of A: 



C(A) = 1.0867- 5.501 x lO^^A 



(9) 



where A is measured in A. Equation |9] was obtained from 
the A-dependent scale-height ratios given in Boselli & 
Gavazzi (1994), by averaging the optically thin and op- 
tically thick cases. In the UV (A ~ 2000 A), this ratio 
is C = 1, so Eq. |8] reduces to a slab model and can be 
numerically inverted: 

t{UV) sec(i) = 0.0259 -f 1.2002 x Ai{UV) + 1.5543 x Ai{UVf 
-0.7409 X Ai{UV)^ + 0.2246 x Ai{UV)'^ (10) 

^'^ Note that when interacting with dust grains, photons can 
be absorbed, scattered out of the line of sight and scattered back 
into it. The term "extinction" refers to the first two processes, 
while "attenuation" encompasses all of them. When talking about 
external galaxies the term "attenuation" is preferred, since it takes 
into account the complex radiative transfer processes resulting form 
the relative geometry of stars and dust. Nevertheless, with this 
caveat in mind we will use both terms interchangeably throughout 
this paper. 
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Once the optical depth in the UV is known, the cor- 
responding value at any other wavelength is given by a 
particular extinction law k{X): 

t(A) = t{UV) X k{X)/k{UV) (11) 

which plugged into Eq. [8] gives us the attenuation at the 
desired band. 

The conversion between the TIR/FUV and TIR/NUV 
ratios into Apuv and AnuV: respectively, is rather in- 
sensitive to the adopted extinction law (Cortese et al. 
2008). Nevertheless, to compute the extinction at other 
wavelengths we must not only choose a particular ex- 
tinction law, but also decide whether to determine A{X) 
by extrapolating from Apuv or Anuv in the equations 
above. Here we use the MW extinction law of Li & 
Draine (2001), assuming Ry = 3.1. Other extinction 
curves are possible, but most of them agree pretty well 
from the near-IR to the NUV bands (Gordon et al. 2003). 
It is beyond the 2175 A bump that large differences arise. 
Therefore, in order to minimize the impact of our par- 
ticular choice of extinction law, we use the NUV band 
rather than the FUV one in Eqs. [TUlandITT] 

5. FITTING PROCEDURE 

In order to find the model that best fits the observed 
multiwavelcngth profiles for each galaxy, a minimiza- 
tion procedure was followed. We generated a grid of mod- 
els with velocities ranging between 130 and 250 km s^^ 
in steps of lOkms"^, plus extra values of 40, 80, 290 
and 360kms~^. As for the spin parameter, we sampled 
the interval 0.02 < A < 0.09 in steps of 0.01; besides, 
we added A = 0.10,0.15 and 0.20 to our grid in order 
to account for possible low surface brightness galaxies 
(LSBs). 

Taking this set of 187 pre-computed models as our 
starting point, we used a 2D interpolation algorithm to 
generate a finer grid of models with steps of 0.001 in A 
and lkms~^ in Vq- We verified that any given property 
of a model galaxy at a certain radius varies smoothly 
enough with A and Vc, so that the corresponding value 
for a model with an intermediate spin and velocity can 
be indeed approximated by means of a 2D interpolation. 

The total of each model was computed by sum- 
ming over data-points at all bands and galactocentric 
distances. By visually inspecting the multi-wavelength 
profiles, we excluded from the fit those radial ranges in 
which the overall emission is dominated by the bulge. 
In those galaxies with sharp outer truncations or anti- 
truncations —which the BPOO models cannot reproduce 
by construction— the outermost regions were excluded as 
wen. This affects NGC 1512, NGC 3621, NGC 4625 and 
NGC 4736, so for these objects the results of the fitting 
only concern the bright inner dislQ . The radial range 
used for the fit in each galaxy is quoted in Table [TJ 

The resulting distribution of values are shown in 
Appendix \K\ If they are to be used to derive confidence 
intervals for the fitted parameters, rather than just to 

More galaxies in the sample exhibit multi-sloped profiles, but 
their breaks are considerably smeared out in the optical and es- 
pecially in the near-IR. Indeed, excluding the outer disks does 
not perceptively change the output of the fit in those cases. For 
NGC 3031, NGC 5194 and NGC 6946, though, we did exclude the 
very outermost parts of the profiles, due to contamination of noisy 
background structures. 



find the best-fitting values, then a proper determination 
of the uncertainties of each data-point and of the models 
is mandatory. In principle, the method assumes that 
any deviation of the observed values with respect to the 
model predictions is entirely due to measurement errors. 
If these errors are properly accounted for when comput- 
ing , then the confidence intervals for the fitted param- 
eters are defined by all models with < Xmin + ^X^j 
where Ax^ depends on the confidence level and the num- 
ber of parameters that are being estimated simultane- 
ously (see e.g. Avni 1976; Press et al. 1992). 

However, we cannot strictly follow this approach in our 
case, because the models do not reproduce the small- 
scale structures of real disks. In an attempt to overcome 
this problem, we first run our fitting code assuming that 
the total uncertainty for each data-point is the quadratic 
sum of the photometric and zero-point errors, plus an 
extra uncertainty of 10%. This additional term serves 
as an initial guess for the intrinsic error of the model, 
and also avoids giving excessive weight to any particular 
band and/or data-point. The typical reduced x^ at this 
stage is of the order of ~ 5. We then compute the relative 
rms of the best-fitting model with respect to the galaxy's 
profiles, both as a function of radius and wavelength. In 
this way we can estimate how well we can expect the 
model to fit that particular galaxy at each band. These 
'error profiles' are then fed to the code in a second run, 
in place of the initial uncertainties. The new reduced 
X^ values are now close to 1, by construction. However, 
the purpose of this two-stage fitting process is not to 
artificially bring the reduced x^ closer to unity, but to 
properly take into account deviations due to small-scale 
features that the models, by construction, are not able 
to reproduce. 

Even after following this process, we found that the 
technique of adding a constant Ax^ offset to the total 
(i.e. not reduced) x^ still yielded unrealistically small 
confidence intervals for Vc and A. A visual inspection 
confirmed that indeed many models outside these confi- 
dence regions were still in very good agreement with the 
observed profiles. Thus, we finally opted for defining the 
boundaries of the confidence intervals with those models 
whose total x^ is twice the minimum one. Therefore, al- 
though the resulting errors in A and Vq cover the range 
of models that visually agree with the galaxy's profiles, 
they are indicative and should not be interpreted in a 
strict statistical way (see also Boselli et al. 2006 in this 
regard) . 

As an example, in Fig. [2] we show the resulting fit for 
the Sc spiral NGC 3198 using the KOI IMF (see the on- 
line edition of the journal for similar plots for the re- 
maining galaxies). The gray data-points show the ob- 
served profiles, corrected only for MW extinction, while 
the black ones are also corrected for internal extinction. 
Both profiles have been deprojected to their face-on val- 
ues by means of the galaxy's morphological axis raticQ. 
The fit is applied to the profiles corrected for internal ex- 
tinction, and only to those points beyond the red dashed 

^■^ Note that, strictly speaking, this deprojection is only valid 
for the profiles corrected for internal extinction, which are the ones 
used in the fit. In the observed profiles, the difference between the 
inclined and face-on values would not just owe to a simple geomet- 
rical projection effect, since the interaction between starlight and 
dust along a different line of sight would also play a role. 
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NGC3198 X=0.063!o.oio V^^igi^g km/s K01 IMF 



NGC2841 A,=0.028!°;oJ8 Vq=271!|§ km/s KOI IMF 
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Fig. 2. — Best-fitting model for the Sc galaxy NGC 3198, using 
the KOI IMF. The gray points show the observed profiles, corrected 
only for Milky way extinction, while the black ones also include 
a correction for the radial variation of internal extinction. Both 
profiles have been deprojected to their face-on values using the 
galaxy's morphological axis ratio. In each panel, the radius along 
the semimajor axis is expressed both in arcseconds (top x axis) and 
in kpc (bottom x axis). The fit is applied to those points in the 
extinction-corrected profiles beyond the dashed red line, in order 
to exclude the bulge. The red curve corresponds to the best-fitting 
model, and the shaded area contains all models with < 2y^- . 



line, which separates the bulge- and disk-dominated re- 
gions of the profiles. In the few cases where we had to 
exclude the outer regions (due to strong up-bendings, for 
instance), the outer limit is marked with a blue dashed 
line. The best-fitting model is shown with a red line, and 
the band with a lighter shade of red contains all models 
with < 2x^i„. 

However, the fit is not always equally good at all wave- 
lengths. In Fig. Owe show the best fitting model for the 
Sb galaxy NGC 2841. Even though the quality of the fit 
is excellent all the way from 4.5 fim to the u band, the 
model overpredicts the luminosity of the galaxy in the 
GALEX bands. This tends to happen mostly in early- 
type spirals, as will be discussed in section [6.2.31 

6. RESULTS 

6.1. Global properties 

The results of the fitting procedure are quoted in Ta- 
ble [2] Prior to going further in our analysis, we will 
briefiy describe the statistical distribution of the model 




Fig. 3.— Same as Fig.[2l but for the Sb spiral NGC 2841. The fit 
is excellent at all wavelengths except in the FUV and NUV bands. 



parameters A and Vq- As mentioned before, for each 
galaxy we have run our fitting code using both the K93 
and KOI IMFs. Given their different content in high- 
mass stars at a fixed total mass, the resulting profiles 
differ in the UV bands, but agree in the optical and near- 
IR ones. The effects of choosing one IMF or another will 
be discussed in detail in section 16.2.31 but for now it 
will suffice to say that neither A nor Vq are significantly 
affected by our particular choice of IMF. Therefore, here- 
after we use the K93 IMF as our default choice, unless 
otherwise mentioned. 

6.1.1. Statistical distribution of the model parameters 

In Fig. |4] we show the resulting histograms of both 
fitting parameters. It can be seen that most galaxies ex- 
hibit values of A and Vc similar to those of the Milky 
Way. In particular, the distribution of rotational veloc- 
ities peaks at 200-220 km s~^. It should be noted, how- 
ever, that neither the SINGS sample nor the smaller sub- 
sample of disks considered here are complete. The well- 
known Schechter (1976) function can be used to fit not 
only the mass and luminosity functions of nearby galaxies 
(Bell et al. 2003), but also their circular velocity distribu- 
tion (Gonzalez et al. 2000; Boissier et al. 2000). In this 
sense, low-mass slow-rotating disks are known to out- 
number more massive and faster-rotating ones. There- 
fore, the velocity distribution shown in Fig. |4] obviously 
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Fig. 4. — Distribution of the circular velocity (a) and spin (b) in our sample. The values adopted for the Milky- Way (A = 0.03 and 
Vc = 220kms~^) have been marked in both panels. The dashed curve corresponds to the probability distribution of A proposed by Mo et 
al. (1998), scaled to match our histogram. 



underestimates the number of low-velocity galaxies that 
would be found in a volume-limited sample. 

Regarding the spin parameter, most disks in our sam- 
ple have A ~ 0.03, the same spin we have adopted for 
the MW. Had we chosen a different MW spin, the result- 
ing distribution of A would have been shifted accordingly. 
Such a peaked histogram is in agreement with a key pre- 
diction of ACDM simulations of galaxy formation: the 
fact that most haloes exhibit the same angular momen- 
tum per unit of mass at any epoch, regardless of their 
total mass and their particular history of mass assembly. 
In this sense, the following analytic expression is known 
to fit the distribution of A obtained in iV-body simu- 
lations simulations (see e.g. Barnes & Efstathiou 1987; 
Warren et al. 1992; Gardner 2001; Bullock et al. 2001; 
Vitvitska et al. 2002): 



p{X)dX 



1 



27rCTA 



exp 



ln2(A/A) 



2al 



dX 

T 



(12) 



In particular. Mo et al. (1998) adopt A = 0.05 and 
a\ = 0.5. Since it is a log- normal function, these val- 
ues should not be understood as the mean and width 
of the distribution. In fact, this function peaks around 
A ~ 0.04, and has a width of ~ 0.05. Interestingly, the 
distribution of spin values of our galaxies agrees well with 
Eq. [121 This implies that even though our sample is not 
representative of a complete one, to some extent it be- 
haves as if it was with regard to A and, by extension, to 
any other quantity that depends primarily on A rather 
than on Vc- 

6.1.2. Comparison with observed values 

Before further proceeding with any detailed analysis, 
we must first verify that the values of A and Vc that 
our fitting code yields for each galaxy agree with the ob- 
served ones. This comparison is not a straightforward 



task in the case of the spin parameter: being a model- 
dependent quantity, it cannot be directly measured in 
real galaxies in the same fashion as the rotational veloc- 
ity. For instance, the scaling laws adopted by Hernandez 
& Cervantes-Sodi (2006) would yield somewhat smaller 
velocities and spins than the BPOO models for galaxies 
more massive than the MW, and viceversa (assuming 
the same reference values of A and Vc for the MW in 
both models). Nevertheless, as commented above. Fig. 2] 
shows that the distribution of A in our sample resembles 
the one usually found in numerical N-body simulations. 

In order to check the accuracy of our circular veloci- 
ties, in Fig. [5^ we compare the theoretical values given 
by the model with the observed rotational velocities re- 
trieved from the Lyon-Meudon Extragalactic DatabaseFl 
(LEDA; Paturel et al. 2003). The latter are determined 
from the width of the 21 cm hydrogen line at different lev- 
els and/or from rotation curves, usually Ha ones. The 
final values provided by LEDA are homogenized and cor- 
rected for inclination. In general, our theoretical val- 
ues for the circular velocity are in agreement with the 
observed ones, which in some sense is expected, given 
that the model incorporates a TuUy-Fisher (TF) relation 
through the adopted ACDM scaling laws (see discussion 
below). 

However, for many galaxies our velocities tend to be 
about 25% larger than those quoted in LEDA. The most 
discrepant outlier in this plot is the Sc spiral NGC 0628, 
for which our fitting code yields Vc = 208kms~^, while 
in LEDA we find a much lower value of 38kms^^. This 
latter velocity is surprisingly small given that NGC 0628 
has an absolute Ks-ha,nd magnitude of —21.64 mags, 
for which one should expect a rotational velocity of 
^ 175 km s^^ according to the TF relation (see below). 
The fact that this galaxy is almost face-on might in- 

1^ |http : //leda ■ unlv-lyonl ■ f r] 
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Fig. 5. — (a): Comparison between the circular velocities derived from the model fitting and the observed ones compiled in the LEDA 
database. Filled symbols correspond to galaxies with a minor to major axis ratio b/a < 0.85, while open ones show almost face-on galaxies 
with b/a > 0.85. (b): Same as panel (a), but using as a reference the circular velocities derived from the THINGS rotation curves. These 
were fitted by Leroy et al. (2008) with an analytic expression, one of whose parameters is the velocity in the flat regime of the rotation 
curve, which is the value we plot here. The small inset shows a comparison between the THINGS velocities and those quoted in LEDA for 
the same objects. 



troduce large uncertainties in the inclination correction, 
thus possibly making the LEDA velocity very uncertain 
for this galaxy. The open circles in Fig.[5K show, however, 
that the LEDA values for many other face-on galaxies in 
the sample agree well with the ones obtained from the 
model. 

There exist in the literature more accurate kinematical 
data for some of our galaxies. Daigle et al. (2006) and 
Dicaire et al. (2008) obtained Ha rotation curves for the 
SINGS galaxies using Fabry-Perot interferometry. Also, 
the SINGS sample overlaps with The HI Nearby Galax- 
ies Survey (THINGS; Walter et al. 2008), for which HI 
rotation curves were derived by de Blok et al. (2008). 
Here we make use of the results of Leroy et al. (2008), 
who parameterized the rotation curves of the THINGS 
galaxies with the following analytical expression: 

Vrot{r) = Wflat [1 - exp (-r/Zflat)] (13) 

where Vrot is the circular rotational velocity at a radius r, 
fflat is the asymptotic velocity where the rotation curve 
is flat, and Zflat is the radial scale length over which Wflat 
is reached. Leroy et al. (2008) derived Wflat and Zfiat from 
the high resolution rotation curves presented in de Blok 
et al. (2008), as well as from the first moment maps for 
those low-inclination galaxies not included in de Blok et 
al. (2008). 

In Fig. we compare the circular velocities derived 
from the model fitting with the Wflat values computed 
by Leroy et al. (2008), for the 17 galaxies we have in 
common. Contrary to what happens with the LEDA ve- 
locities in Fig. [S^, no systematic shift appears. Our ve- 
locities are within 20-25% from the THINGS ones. What 
is more, this scatter is uncorrelated with the one between 
the predicted and observed photometric profiles. The ve- 



locities quoted in LEDA for some galaxies in the sample 
of Leroy et al. (2008) tend to be lower than the values de- 
rived by those authors, and part of the discrepancy seems 
to be associated with differences in the corrections for 
inclination. Given the exquisite quality of the THINGS 
data and the homogeneity in the derivation of the ro- 
tation curves, we conclude that most of the systematic 
offset seen in Fig. [5^ is likely an issue of the LEDA values. 

Another way to check the validity of our model rota- 
tional velocities consists of trying to reproduce the TuUy- 
Fisher relation (TuUy & Fisher 1977). This tight empiri- 
cal relation links the intrinsic luminosity of a galaxy with 
the amplitude of its rotation curve. The former quantity 
traces the stellar mass, while the latter probes the total 
gravitational mass. Therefore, any successful model of 
disk evolution must be able to reproduce this observed 
correlation. From a observational point of view, numer- 
ous studies have shown how the slope and zero-point of 
the TF relation vary with wavelength in the optical and 
near-IR, mainly due to changes in the mass-to-light ratio, 
together with extinction if it is not properly accounted 
for (Bell & de Jong 2001; Courteau et al. 2007; Pizagno 
et al. 2007; Blanton & Moustakas 2009). 

The scaling laws adopted by BPOO imprint a built- 
in TF relation in the models through Eq. However, 
once star formation is implemented in a self-consistent 
way, the resulting slope and zero-point of the TF rela- 
tion might vary with wavelength due to changes in the 
mass-to-light ratio (see also Ferreras & Silk 2001). When 
comparing the TF relation resulting from the models 
with several empirical ones in the / band from different 
authors, BPOO found a good agreement, although their 
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Fig. 6. — Absolute ii'g-band magnitude of our galaxies as a func- 
tion of the circular velocity of the corresponding best- fitting model. 
The absolute magnitudes have been derived from the asymptotic 
values, and are corrected for internal extinction. The shaded bands 
corresponds to the empirical TF relations of Verheijen (2001) and 
Courteau et al. (2007). The width of the bands shows the observed 
1(T scatter of those fits. The solid and dashed lines show direct and 
reverse linear fits to our data, respectively. 



theoretical TF relation yielded somewhat larger veloci- 
ties for a given absolute /-band magnitude. 

In Fig. ini we plot the iCs-band absolute magnitude of 
our galaxies as a function of the circular velocity re- 
sulting from the model fitting. The absolute magni- 
tudes were computed from the asymptotic values pre- 
sented in Paper I, and were corrected for internal extinc- 
tion from the global TIR/UV ratio as described above. 
The resulting median extinction in the Ks band is just 
~ 0.03 mag. We compare our results with the empiri- 
cal Ks-hand TF fits of Verheijen (2001) and Courteau 
et al. (2007). The model TF relation lies slightly below 
the observed ones, with values of Vq roughly 0.05-0.1 dex 
larger for a given absolute magnitude, especially for the 
most massive galaxies. This translates into a relative off- 
set of 10-20%, which is much lower than the systematic 
scatter seen in Fig. [SK- After applying a direct and a re- 
verse weighted linear fit to our data-points, the following 
relations are obtained: 

Mks 1.3 ± 1.4 - (10.01 ± 0.60) X log Vq (14) 
(rms — 0.61 mag) 
log Vc = -0.04 ± 0.10 - (0.1069 ± 0.0049) x Mk115) 
{rms — 0.068 dex) 

where is expressed in AB magnitudes and Vq in 

kms~^. Here we rely on the Kg band rather than on 
optical ones in order to minimize the effects of internal 
extinction, but comparisons between models and obser- 
vations at other wavelengths can be found in BPOO. 

6.1.3. Trends along the Hubble sequence 

It is illustrative to discuss whether the derived values 
of A and Vq depend on the morphological type. In Fig. [7] 




Fig. 7. — Dependence of the spin parameter (a) and the circu- 
lar velocity (b) on the morphological type T. Given that the T 
types are integer values, small random horizontal offsets have been 
applied for the sake of clarity. 



we plot the model parameters as a function of the Hubble 
type. The rotational velocity is clearly correlated with 
the morphology of the galaxies, with early-type disks 
rotating faster —and hence being more massive— than 
late-type ones (Bosma 1978; Roberts 1978; Rubin et al. 
1985). However, there is no apparent trend between the 
spin parameter and the Hubble type, except maybe an 
increased scatter in Sdm-Sm galaxies. 

Although the fact that the Hubble type does not seem 
to depend on the spin parameter might look surprising 
at first glance, it supports the findings of numerical sim- 
ulations. The number density of galaxies per velocity 
interval can be parameterized with a generalized form of 
the Schechter function, resulting from the combination 
of the luminosity function and the TF relation at a given 
band (Gonzalez et al. 2000). This distribution is obvi- 
ously monotonic, with massive disks being rather scarce 
compared to low-mass ones. Given that Vq depends on 
the Hubble type, the probability distribution of A would 
be also monotonic if this parameter was also correlated 
with the morphological type. However, we have seen that 
numerical N-body simulations conclude that most haloes 
usually exhibit the same 'universal' spin value A 0.04 
quite irrespective of their mass. Therefore, the lack of 
correlation between A and the Hubble type agrees with 
this result. 

6.2. Potential sources of systematic errors and 
discrepancies in the UV 

In this section we will explore different issues that could 
potentially introduce biases in our results. We will pay 
special attention to those mechanisms that could be re- 
sponsible of the excess in the UV luminosity predicted 
by the model in early- type disks (see Fig. [3]). 
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6.2.1. Inclination 

After the profiles have been corrected for internal ex- 
tinction, they are also deprojected prior to fitting them 
with the models. Correcting for inclination modifies the 
overall surface brightness level, so the resulting values 
of A and Vq could be potentially affected by the precise 
value of the inclination angle. For consistency with the 
way in which the surface photometry was done, the de- 
projection was done by means of the observed axial ratio, 
q — b/a^ so that cosi — q. However, if disks are mod- 
eled as oblate spheroids, then the inclination angle can 
be corrected for the effects of the intrinsic thickness of 
the disk: 

COs2z' = (g2_q2)(^_^2)-l (jg) 

Here go is the intrinsic axial ratio of the disk when ob- 
served edge-on. It is customary to assume a constant 
value of go — 0.2 (see e.g. Courteau et al. 2007), although 
some authors argue that this value could decrease to- 
wards go — 0.1 in late- type spirals (Haynes & Giovanelli 
1984; Dale et al. 1997). Taking into account the thickness 
of the disk slightly increases the inclination angle. As a 
result, the surface brightness /i' deprojected via Eq. ll6l is 
fainter than the value /i deprojected assuming cosi ~ q: 



/ , cos i , / 1 

Au = u' - u = -2.5 log = -2.5 log J — 

COSI y 

'(17) 

It can be seen that A/i lies well below 0.01 mag for 
i < 50°. At the largest inclination angles in our sample 
(i ~ 80° for a couple of objects) the offset in the surface 
brightness reaches 0.2 mag. By inspecting the output 
profiles of the BPOO models, we determined that A/i ~ 
0.2mag corresponds to a decrease of just lOkms"^ in 
Vc with respect to the case in which the disk thickness 
is ignored. The spin parameter may increase slightly to 
compensate for the decrease in the radial scale-length 
due to the lower circular velocity (see Fig. [T|). 

The most inclined spirals in our sample are NGC 0024 
(Sc; b/a = 0.22) and NGC 7331 (Sb; b/a ^ 0.35). Our 
fitting code yields Vc = 93^1 knis-^ and A = 0.07lj:^:EJi^ 
for NGC 0024, if we assume that go = 0.13 (Dale et 
al. 1997), whereas we get Vc = 98^5 km s"^ and A = 
0.067;^ |^i^ when go = 0. As for NGC 7331, we obtain 
Vc = 250ti^kms-i and A = 0.064^°;°^^ when assuming 
go = 0.20 (Dale et al. 1997), whereas ignoring the thick- 
ness leads to Vc = 263^16 km s^^ ^nd A = 0-059toml- 
Therefore, employing a non-zero value of go in disks close 
to edge-on leads to differences of merely ~ 5% in both 
A and Vc, which lie within the estimated uncertainties. 
The vast majority of the galaxies in our sample are not 
so inclined as the two extreme examples discussed here, 
so including or not the disk thickness has an entirely 
negligible impact in the model parameters. 

6.2.2. Dust attenuation 

Inaccuracies in the internal extinction correction might 
be responsible of the large offset between the observed 
and predicted UV profiles that we find in some disks. As 
explained in Sectional we first compute Apuv and j4nuv 
from the TIR/FUV and TIR/NUV ratios, respectively. 
The extinction in the optical and near-IR bands is then 



derived from Anuv assuming a sandwhich geometry and 
a MW-like extinction curve. 

By means of radiative transfer models, several studies 
(Buat & Xu 1996; Meurer et al. 1999; Gordon et al. 2000; 
Witt & Gordon 2000; Buat et al. 2005) have shown that 
the TIR/UV ratio constitutes a robust proxy for the in- 
ternal extinction in galaxies. In particular, it appears to 
be quite insensitive to the relative distribution of stars 
and dust and to the extinction curve (see, e.g., Fig. 12b 
in Witt & Gordon 2000). In the particular case of the 
recipes of Cortese et al. (2008) —the ones employed here 
to compute Apuv and Anuv~: variations in the dust 
geometry and the extinction law modify the derived at- 
tenuations by less than 0.2 mag (see their Fig. 9). 

When computing the extinction A\ at other wave- 
lengths by extrapolating from Anuv, our particular 
choice of dust geometry end extinction curve will indeed 
play a role. However, in the case of the near-IR bands 
the impact will be minimal: at those wavelengths the 
internal extinction will be still close to zero even if the 
sandwhich model is not appropriate or if the galaxy's 
extinction curve does not resemble that of the MW. In 
brief, for a given galaxy where the model that best fits 
the near-IR profiles significantly overestimates the UV 
ones, switching to another dust-to-stars geometry or ex- 
tinction law would not fix the problem. The extinction in 
the UV would not vary noticeably due to the robustness 
of the TIR/UV ratio against these changes, and the ex- 
tinction in the near-IR would remain close to zero. Vari- 
ations would show up at intermediate wavelengths, in 
the optical range, but the model seems to match the op- 
tical profiles even in those cases where it fails in the UV 
(Fig.©. 

Another aspect of the internal extinction correction 
that should be taken into account is the relative role of 
young and old stars in heating the dust. In galaxies with 
low SFR per unit of mass (specific SFR, or sSFR), the 
contribution of old stars to the dust heating will be larger 
than in more actively star-forming galaxies, thus bias- 
ing the energy balance argument on which the TIR/UV 
method relies. Therefore, a given calibration between 
TIR/UV and A\jy that is valid for late-type spirals will 
overestimate the actual extinction in early-type ones. 

In Paper II we used two different recipes to derive at- 
tenuation profiles for the SINGS galaxies: the calibration 
of Buat et al. (2005), which does not account for the 
varying extra dust heating due to old stars, and the age- 
dependent calibration of Cortese et al. (2008), the one 
finally adopted here. The difference between both pre- 
scriptions is negligible in disk-dominated galaxies. How- 
ever, in SO/a-Sab galaxies and the bulges of later- type 
spirals, the method of Buat et al. (2005) yields extinc- 
tion values ~ 0.5 mag larger than those obtained with the 
recipes of Cortese et al. (2008), since the extra contribu- 
tion of old stars to the dust heating is not "removed" 
in the former. Therefore, using the calibration of Buat 
et al. (2005) would bring the observed profiles closer to 
the model predictions. However, this better agreement 
would be achieved at the expense of using an extinc- 
tion recipe that is known not to be valid in early-type 
spirals. Nevertheless, the difference of ~ 0.5 mag is not 
large enough to account for the UV discrepancies found 
in many of our disks. 
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Fig. 8.— (a): The IMFs of Salpeter (1955), Kroupa et al. (1993) 
and Kroupa (2001), normalized to a total stellar mass of IMq. 
(b): Ratio of the K93 and KOI IMFs to the Salpeter one. 



6.2.3. Initial Mass Function 

The IMF, usually denoted as $(M), indicates the num- 
ber of stars in the mass interval M to A/ -f dM formed 
in a given burst. The classical IMF of Salpeter (1955) 
consists of a single power-law across the whole stellar 
mass range, $(M) oc m~", where a = 2.35. In this work 
we have used two grid of models with the K93 and KOI 
IMFs, which are multi-sloped: 



0.1 < M/Mq < 0.5 
0.5 < M/Mq < 1.0 
1.0 < M/Mq 



«K93 = 1-3 ; axoi = 1-3 
«K93 = 2.2 ; ttKoi = 2.3 
aK93 =2.7; aKoi -2.3(18) 



These IMFs are shown in Fig. |8l where they have been 
normalized such that: 



<S>{M)MdM = 1 



(19) 



where the integration is carried out between Mmin ~ 
0.1 M© and M^ax = lOOMgE!!. The KOI IMF has more 
short-lived massive stars for a given total stellar mass, so 
it yields larger UV fluxes for the same near-IR and opti- 
cal luminosities. Therefore, choosing one IMF or another 
obviously affects the ability of the model to simultane- 
ously fit all the multi- wavelength profiles. 

In order to quantify the discrepancy between the model 
predictions and the actual profiles, for each galaxy and 
wavelength we have computed the average difference be- 
tween the observed and model surface brightness within 

^'^ Note that while the Salpeter IMF is normally assumed to hold 
over the range 0.1-100 Mq, the original Salpeter (1955) study only 
included stars in the range 0.4-10 Mq 
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Fig. 9. — Average difference between the observed and model 
surface brightness profiles for each galaxy at different wavelengths. 
The average offset was computed within the radial range used dur- 
ing the fitting procedure. Note that the "observed" profiles are 
already corrected for internal extinction. 



the radial range considered in the fitting procedure. "Ob- 
served" is used here in contrast to "model" , but the "ob- 
served" profiles we consider are of course those corrected 
from internal extinction. The results are shown in Fig. [S] 
for each IMF and morphological type. 

It is clear that neither IMF provides a satisfactory fit 
at all wavelengths and for all morphological types simul- 
taneously. The largest discrepancies are found in the 
UV range, as expected. The KOI IMF yields an excel- 
lent fit for most Sc-Sd spirals, as well as for some Sb-Sbc 
ones. However, it significantly overestimates the UV fiux 
in early-type disks. The K93 IMF, on the other hand, 
partly mitigates the problem in some Sb-Sbc galaxies, 
but the perfect agreement between model and observa- 
tions found in Sc-Sd disks is partially lost. Very late-type 
disks, however, are better fitted with the K93 IMF than 
with the KOI one. 

Elucidating whether a varying IMF is the actual rea- 
son behind these discrepancies is not straightforward, as 
other mechanisms could introduce similar systematic bi- 
ases in the UV profiles. For instance, the UV bands 
can be affected by variations in the SFR taking place on 
timescales of ~100 Myr, whereas the u band will respond 
to variations over ^500 Myr. This increase in the char- 
acteristic timescale between the UV and optical bands 
seems to be indeed rather sharp (Boissier et al. 2011, in 
prep) . Therefore, short-scale SFR variations could partly 
explain why the discrepancies between the observed and 
predicted profiles drops abruptly between the UV and u 
bands. Also, the calibration of the stellar atmospheres 
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Fig. 10. — Comparison of the values of the spin parameter and circular velocities obtained when using the IMFs of Kroupa et al. (1993) 
and Kroupa (2001). The solid line shows the corresponding 1:1 relations. 



used in the spectral synthesis models might not be as 
accurate in the UV as they are in the optical, since ob- 
servations of stars in the UV are not so easily accesible. 

It has been argued that the so-called Integrated Galac- 
tic Initial Mass Function (IGIMF), which takes into 
account the clustered nature of star formation, varies 
among galaxies, leading to a scarcity of massive stars in 
galaxies with low global SFRs (Weidner & Kroupa 2005). 
Krumholz & McKee (2008) brought forward a physical 
explanation to the apparent dearth of massive stars in 
low-density regions, arguing that gas column densities of 
at least lgcm~^ are required to halt cloud fragmenta- 
tion and produce massive stars. Stochastic formation of 
high mass stars might also play a role in galaxies with 
very low levels of star formation (Lee et al. 2009). 

These mechanisms have been proposed to explain the 
observed discrepancies between the UV and Ha-derived 
SFRs in dwarf galaxies (Pflamm-Altenburg et al. 2007, 
2009; Lee et al. 2009) and in low surface brightness galax- 
ies (Meurer et al. 2009). In this regard, even though both 
the K93 and KOI IMFs are canonical, it would make 
sense that Sdm-Sm disks are better fitted with a top- 
light IMF like the K93 one rather than with the KOI 
one. However, the effects described here seem to arise 
mainly in irregular galaxies with much lower levels of SF 
than the disks considered in this work (see also Boselli et 
al. 2009). In any case, we cannot appeal to these mecha- 
nisms to explain the discrepancies in the UV bands found 
in early- type disks. 

Nevertheless, it is worth asking whether our particular 
choice of IMF affects the results of the fitting. Thanks 



to our two-step fitting procedure, which estimates the 
intrinsic error of the model for each galaxy and band af- 
ter a first run, the UV bands are automatically assigned 
a relatively large error whenever large discrepancies are 
found. This prevents the UV bands from biasing the fit 
at longer wavelengths. In Fig. [10] we compare the values 
of A and Vc obtained with both IMFs. The differences 
are obviously negligible and much smaller than the es- 
timated uncertainties, since the fit is still excellent in 
the optical and near-IR bands even when it fails in the 
UV. Therefore, even if the very recent level of SF (either 
massive or total) predicted by the models for the last 
few hundreds of Myrs is not entirely reliable, the overall 
star formation history across longer timescales can still 
be trusted. 

Besides modifying the emitted UV fiux, changing the 
IMF also has a significant impact in the resulting metal- 
licity profiles predicted by the model. In Fig. [11] we 
compare the central abundances and radial gradients ob- 
tained with both IMFs. Given that the KOI IMF is richer 
in massive stars than the K93 one, it produces more Type 
II supernovae and metals, leading to oxygen abundances 
which are ~ 0.62 dex larger than those resulting from the 
K93 IMF. Note that the uncertainties associated with 
different metallicity calibrations are typically ~ 0.3 dex 
(see Moustakas et al. 2010). The radial gradients, on 
the other hand, remain nearly unchanged, the KOI ones 
being just mildly steeper. 

We must also check whether the oxygen abundance 
profiles predicted by the model are in agreement with 
the observed ones. As in Paper II, here we rely on the 
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Fig. 11. — Comparison of the central abundances (a) and gradients (b) obtained when using the IMFs of Kroupa et al. (1993) and Kroupa 
(2001). The sohd hnes in each panel correspond to the 1:1 relation, while the dashed lines show linear fits where only the zero-point has 
been left as a free parameter. The corresponding offset is quoted in each panel. 



metallicity zero-points and gradients measured by Mous- 
takas et al. (2010), using the calibration of Kobulnicky 
& Kewley (2004). In Fig. [12] we plot the central oxy- 
gen abundances as a function of the predicted values 
using both IMFs. The values computed with the K93 
IMF are just 0.1 dex larger than the observed central 
abundances, whereas those yielded by the KOI IMF are 
roughly 0.7 dex larger. Such high metallicites have never 
been observed; moreover, the flat slope for massive stars 
of the KOI IMF would lead to a significant depletion of 
deuterium through astration in the solar neighborhood. 

As pointed out by several authors (see, e.g., Kewley 
& Ellison 2008; Moustakas et al. 2010), different meth- 
ods used to measure the gas-phase metallicity from ob- 
served emission-line spectra can yield largely discrepant 
results. Besides the calibration of Kobulnicky & Kewley 

(2004) , Moustakas et al. (2010) also derived metallicity 
profiles using the calibration of Pilyugin & Thuan (2005). 
Both methods are based on the strong-line R23 parame- 
ter (Pagel et al. 1979). The empirical method of Pilyu- 
gin & Thuan (2005) is calibrated on HII regions having 
direct abundance measurements based on the electron 
temperature. The theoretical calibration of Kobulnicky 
& Kewley (2004), on the contrary, relies on photoioniza- 
tion models. The method of Pilyugin & Thuan (2005) 
leads to oxygen abundances '^O.Bdex lower than those 
obtained with the Kobulnicky & Kewley (2004) calibra- 
tion. Therefore, had we used the Pilyugin & Thuan 

(2005) values in Fig. [T^l even the metallicites predicted 
with the K93 IMF would be too high compared with the 



observed ones, and the discrepancy with the KOI IMF 
would be even larger than it already is. As thoroughly 
discussed by Moustakas et al. (2010), neither empirical 
nor theoretical strong-line methods are devoid of prob- 
lems. On one hand, empirical methods like the Pilyugin 
& Thuan (2005) one might fail outside the metallicity 
range spanned by the HII regions used in the calibrations, 
and can therefore underestimate the true abundances in 
the high-metallicity regime. On the other hand, theoret- 
ical methods such as the Kobulnicky & Kewley (2004) 
one adopt several simplifying assumptions on the prop- 
erties of the gas clouds and the ionizing stars that can 
make the derived metallicities higher than the actual 
ones. The two calibrations discussed here bracket the 
range of abundances usually obtained with other meth- 
ods, so the actual metallicities probably lie somewhere 
in between. 

As for the radial gradients, they are plotted in Fig.fTSl 
both in units of dexkpc"^ (upper panels) and normalized 
by the optical radius R25 (lower panels). The gradients 
predicted by the model are always steeper than the ob- 
served ones, although for most galaxies the difference is 
just about ^ 0.015 dex kpc^^. When expressed in terms 
of the optical radius, the model predicts a roughly con- 
stant gradient of -0.8dex/R25 for the K93 IMF, and 
-0.9dex/R25 for the KOI one. This "universal" oxy- 
gen gradient, which does not depend on galaxy mass, is 
usually found in observations (see e.g. Henry & Worthey 
1999). However, the spread is larger than that predicted 
by the models, with observed gradients ranging from 
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Fig. 12. — Comparison between the observed and predicted central abundances for those SINGS galaxies with metallicity profiles computed 
by MoustaJjas et al. (2010). The average offsets for the K93 and KOI IMFs are shown with a dashed and dot-dashed line, respectively, 
while the solid one corresponds to the line of equality. 
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Fig. 13. — Comparison between the observed and predicted abundance gradients for those galaxies in the sample of Moustakas ct al. 
(2010). Panels (a) and (b) show the gradients in units of dexkpc"^ for both IMFs, whereas in panels (c) and (d) the gradients are expressed 
in units of the optical radius R25. Galaxies are sorted out into non-barred (A), barred (B) and mixed (AB). The solid line marks the 1:1 
relation in all cases. 
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— ldex/R25 to almost flat gradients. The fact that the 
comparison of gradients in units of dexkpc"^ is some- 
what better than in dex/R25 is not surprising, since the 
former are mainly determined by the radius of the disks, 
which varies substantially within the sample. Therefore, 
the comparison between observed and model gradients 
in dexkpc"^ is more a test of the prediction of disk sizes 
than of abundance gradients. 

Bars may induce radial gas flows that can yield shal- 
lower metallicity gradients (see, e.g., Martin & Roy 
1994). Two out of the three barred galaxies in Fig. [13] 
are among the ones with the flattest observed gradients, 
although the subsample considered here is obviously too 
small to extract any statistically significant conclusion. 
Recent numerical N-body simulations have shown that 
radial stellar migration can also fiatten the final metal- 
licity profiles (Roskar et al. 2008; Martinez-Serrano et al. 
2009; Sanchez-Blazquez et al. 2009), although these re- 
sults refer to the stellar metallicity and not the gas-phase 
one. 

To summarize this dicussion on the IMF, most Sc-Sd 
spirals require the KOI IMF in order not to underesti- 
mate the UV luminosity, but at the expense of ending 
up with oxygen abundances much larger than the ob- 
served ones. Conversely, the K93 IMF provides a some- 
what better fit for the remaining Hubble types, and is 
able to reproduce the correct present-day abundances. 
Whether the better agreement with the observed metal- 
licities should be given a special importance is unclear, 
since neither the measured values nor the predicted ones 
are devoid of possible sources of large systematic errors. 
On one hand, the different existing calibrations used to 
compute the oxygen abundance from observed line ra- 
tios may lead to systematically different values, as ex- 
plained before (Kewley & Ellison 2008; Moustakas et al. 
2010). On the other hand, uncertainties in the stellar 
yields used in the disk evolution models will affect the 
predicted metallicities. Besides, we cannot neglect the 
possibility that the library of synthetic spectra used to 
compute the multiwavelength profiles could be quite off 
in the UV range. 

With these issues in mind, we opt for a compromise 
solution and adopt the K93 IMF as our default choice in 
this paper, otherwise stated. This choice is done without 
any prejudice to the KOI which, we reiterate, provides 
an excellent fit for many disks in our sample, especially 
Sc-Sd ones. 

6.2.4. Gas content and SF efficiency 

Apart from the broadband photometric profiles, the 
model also outputs the total gas density profiles, which 
can be therefore used as another observational constrain 
of the accuracy of the model. In particular, and con- 
tinuing the discussion presented in the previous section, 
the gas fraction might shed some light into the disagree- 
ment in abundance gradients. The baryonic gas frac- 
tion is defined as fg — AIg/{AI^, + Mg), where and 
Mg are the stellar and gas masses, respectively. As star 
formation proceeds, gas is progressively consumed and 
transformed into stars, which later enrich the remaining 
gas with heay elements. Therefore, fg is expected to de- 
crease with 0/H, although the particular trend might 
depend on the presence of inflows or outfiows (Boissier 
et al. 2001; Garnett 2002). 




8 8.5 9 9.5 10 

12+log(0/H) 

Fig. 14. — Trend between the gas fraction and the gas-phase 
oxygen abundance. The shaded bands show the model predictions 
for the K93 and KOI IMFs. The black lines show the observed 
gas fractions of our galaxies, resulting from the combination of 
atomic and molecular gas profiles with stellar mass ones. The latter 
were derived from the observed 3.6 /^m profiles, adopting the M/L 
ratio yielded by the model at this band. Solid lines indicate those 
galaxies with H2 profiles, whereas dashed ones are used when only 
HI profiles were available. The observed gas profiles were multiplied 
by 1.36 to include helium and heavy elements. The model profiles 
already show the total gas surface density. 

In Fig. [13] we compare the gas fraction profiles pre- 
dicted by the model with the actual ones found in our 
galaxies. According to the model, fg decreases with Vc 
and increases with A (Boissier et al. 2001). The shaded 
areas shown here cover the range of values of Vc and 
A found in our sample. As for the observed quantities, 
the atomic gas profiles were measured on HI maps from 
THINGS, using the same ellipticity, position angle and 
radial step as those used when computing the UV, optical 
and near-IR profiles. Whenever possible, we included CO 
profiles compiled from the literature to account for the 
molecular gas content. The CO profiles were converted 
into H2 ones by means of a metallicity-dependent Xco 
conversion factor (Boselli et al. 2002), using the local 
oxygen abundance at each galactocentric distance. Solid 
lines are used for those galaxies with available molecular 
gas profiles (although most if not all of the gas content 
in the outer parts is entirely atomic). Dashed lines corre- 
spond to galaxies for which only HI data were available. 
In all cases, we multiplied the hydrogen profiles by 1.36 to 
account for helium and heavy elements. Stellar mass pro- 
files for each galaxy were derived from the 3.6 /um ones, 
using the M/ L ratio yielded by the model, and combined 
with the gas profiles to obtain fg{r). 

As already discussed before, the KOI IMF yields ex- 
cessively large metallicities, whereas the K93 provides a 
much better agreement with the observed values. While 
the model is successful at reproducing the observed 
trend, for some galaxies fg seems to decrease faster with 
0/H compared to the model predictions, consistently 
with the observed 0/H gradients being fiatter than the 
model ones. 

The gas profiles can also shed light on the discrepancy 
between the observed and predicted UV profiles, which 
may arise from differences in the amount of gas that set- 
tles onto the disk or the way it is transformed into stars. 
In Fig. [15] we show the ratio between the gas surface 
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Fig. 15. — Ratio between the observed gas surface density and 
the value predicted by the model at different radii for each galaxy 
(see text for details). Solid lines correspond to galaxies with both 
atomic and molecular gas profiles available, whereas dashed ones 
are used for those objects laking CO profiles. The observed gas 
profiles have been multiplied by 1.36 to account for helium and 
heavy elements. The model profiles already include the total gas 
surface density. 



density predicted by the model and the actual one for 
each galaxy as a function of radius. Here we plot the 
model gas density obtained with the K93 IMF; the val- 
ues yielded by the KOI IMF are within ±20% of the ones 
shown here. 

In the disks of Sc-Sd spirals the gas profiles predicted 
by the models agree with the observed ones. However, in 
general the model overestimates the actual gas content 
by a factor of < 10, which is considerable. Does this 
lead to systematic differences in the emerging FUV lu- 
minosity? In Fig. [in] we compare the offset between the 
observed and predicted FUV brightness with the corre- 
sponding observed-to-predicted gas ratio at each radius. 
While there is a non- negligible scatter, a clear trend is 
seen, in the sense that the model overestimates the FUV 
brightness wherever it also predicts too much gas con- 
tent. 

It is worth investigating whether the star formation 
efficiency plays a role in driving the predicted UV pro- 
files away from the actual ones. After all, even if the 
gas profiles are correctly reproduced by the model, the 
emitted UV light will still depend on the particular way 
in which the model handles the process of star forma- 
tion. As explicitly stated in Eq. [1] the model relies on 
a hybrid star formation law that combines the classical 
Schmidt-Kennicutt law with an orbital term (see Leroy 
et al. (2008) for an extensive analysis of various other 
laws) . The value of the efficiency a in Eq. [1] was cho- 
sen to reproduce the present-day observed gas fraction 
in nearby galaxies (Boissier et al. 2003). 

Even though the efhciency a is kept constant in the 
modelj*^. it may actually vary among galaxies or even 
within them. This assumption can be tested, given that 
a depends on the SFR density, the gas density and the 

Note that although a can be interpreted as the efficiency of 
the SFR, by star formation efficiency (SFE) one usually means 
the ratio of the SFR and the gas mass. In our case, SFE{t, r) = 
Q:Sg(t, r)'''^y(r)r~^, so even under the assumption of a constant 
value of a, SFE is still a varying function of radius and time. 
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Fig. 16. — The observed-to-predicted gas ratio is plotted against 
the discrepancy between the actual UV surface brightness —already 
corrected for internal extinction— and the value predicted by the 
model. This UV discrepancy is expressed in terms of a magnitude 
offset (right y axis) and the corresponding flux ratio (left y axis). 
Unlike Fig. (9] where this offset was averaged all over the disk, here 
we plot the full range of values at different galactocentric distances, 
so each galaxy follows a track in this plot. 



rotational velocity at each radius, all of which can be 
measured. In Fig. [T7] we plot the radial variation of a 
for each galaxy. The SFR profiles were computed from 
the extinction-corrected FUV ones, using the conversion 
factor of Kennicutt (1998). This recipe assumes solar 
metallicity, a constant SFR over the last few Myrs and a 
Salpeter (1955) IMF. The small tics at the bottom right 
show how the profiles whould shift with the K93 and KOI 
IMFs. The gas surface density was again obtained from 
THINGS and CO profiles. Finally, the circular velocity 
at each radius was derived from the fits of Leroy et al. 
(2008) to the THINGS rotation curves (Eq. For 
the few THINGS galaxies not considered by Leroy et al. 
(2008), the rotation curve from the best fitting model 
was used instead as a proxy for the actual one. 

Fig. [T71 shows that, except for a few objects, a remains 
roughly constant with radius within most disks, although 
it does vary from galaxy to galaxy. The empirical val- 
ues of a are consistent with or somewhat larger than the 
constant value assumed in the model, shown here as a 
horizontal line. The average scatter is nonetheless con- 
siderable, around 1 dex. There are several factors con- 
tributing to this dispersion in the measured values of 
a. On one hand, variations in the IMF will affect the 
adopted calibration for the SFR as a function of the FUV 
luminosity. Moreover, the assumption of a constant SFR 
over the last few Myrs might not hold in some cases. 
Besides, departures from solar metallicity —both among 
and within galaxies— will also affect the SFR calibra- 
tion. Uncertainties in the CO-to-H2 conversion factor 
will modify the total gas surface density, although this 
is only a concern in the innermost regions. Neverthe- 
less, part of the observed scatter likely reflects intrinsic 
variations of a among galaxies. 

Interestingly, when plotting the observed-to-predicted 
UV offset as a function of the empirical values of a we 
do not observe any significant trend between both pa- 
rameters, in contrast with Fig. I16[ where the data cloud 
was clearly tilted. In other words, the UV discrepancy 
appears to be independent of the precise value of a, so a 
morphology-dependent value of a cannot account for the 
mismatch between the observed and theoretical UV pro- 
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Fig. 17. — Empirical value of the efficiency parameter a in Eq. [l] 
(see text for details on how a. was derived from observations). The 
units of a are such that Ssfr('') is measured in Mq pc~^ Gyr~^, 
SgasCf") in Mq pc~^, r in kpc and V{r) in kms"'^. The horizontal 
solid line marks the constant value assumed in the models (see 
Boissier et al. 2003). SsFR(r) was derived from the extinction- 
corrected UV profiles using the recipe of Kennicutt (1998), which 
assumes a Salpeter (1955) IMF. The small tics at the bottom right 
show how the plotted curves would shift with the K93 and KOI 
IMFs. 

files. Therefore, the bottom line of this analysis is that 
the discrepancies in the UV profiles of early-type disks 
are mainly due to an excessive amount of gas retained by 
the model in these galaxies, rather than to its subsequent 
conversion into stars. 

Some of the galaxies in our sample with the largest UV 
discrepancies belong to the Virgo cluster or the Coma 
Cloud. Hydrodynamical interactions with the hot inter- 
galactic medium or gravitational ones with other cluster 
members might have removed part of the gas in their 
disks, thus quenching the recent star formation activity 
(see Boselli & Gavazzi 2006 for a review). This effect 
would lead to significantly lower UV fluxes than those 
predicted by the model. In order to test whether this is 
actually the case, we have computed the so-called HI de- 
ficiency for the galaxies in our sample. This parameter, 
defined by Haynes & Giovanelli (1984), compares the ob- 
served HI mass of a galaxy with the typical HI mass of 
isolated field galaxies with a similar morphological type 
jnobs linear optical diameter I?°pt: 

HI def = (log MHI(T°•^^ Bf^C)) " log Mm^ (20) 

Therefore, positive HI deficiencies correspond to galaxies 
with lower gas contents than those of similar but isolated 
objects. Here we use the calibration of MHi(r°''% -D°pt) 
derived by Solanes et al. (1996), since it relies on a larger 
sample of galaxies than that used in the seminal paper 
by Haynes & Giovanelli (1984). Following Solanes et al. 
(2001), the reference HI masses for galaxies later than 
Sc have been computed following the prescriptions for 
Sc ones. The intrinsic scatter in measurements of the HI 
deficiency is typically ±0.2-0.3 dex (Haynes & Giovanelli 
1984; Solanes et al. 1996), and it is customary to assume 
that galaxies with deficiencies lower than 0.03-0.05 posses 
normal HI contents. 

In Fig. [18] we plot the average UV offset between the 
observed and predicted UV profiles as a function of the 
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Fig. 18. — Offset between the observed and model UV surface 
brightness as a function of the HI deficiency. The latter is obtained 
by comparing, in logarithmic scale, the total HI mass of a given 
galaxy and the typical HI mass of isolated galax ies with similar 
morphological types and linear diameters (see Eq. I20II . 



HI deficiency. The integrated HI masses of the SINGS 
galaxies have been taken from Draine et al. (2007). There 
are indeed some galaxies with simultaneously large HI 
deficiencies and UV offsets. In NGC 4569, NGC 4579 
and NGC 4826, the model perfectly fits the optical and 
near-IR profiles, but it substantially overestimates the 
UV ones, especially in the outer regions. 

NGC 4826 is an anemic spiral in the Coma 1 Cloud 
(van den Bergh 1976; Boselh & Gavazzi 2009). It is 
known to host two counter-rotating gaseous disks that 
possibly point toward a past merger event (Braun et al. 
1992). It also exhibits a central dark lane that has earned 
this object the nickname "The black-eye galaxy". Its 
stellar component extends beyond r > 13 kpc, as shown 
by the optical and near-IR profiles, but most of the star 
formation activity is currently restricted to r < 3 kpc, 
where most of the gas is located. 

Something similar happens with NGC 4569 and 
NGC 4579, two other anemic spirals in the Virgo Cluster 
(van den Bergh 1976). In order to properly fit the multi- 
wavelength profiles of NGC 4569, Bosehi et al. (2006) 
employed a modified version of the BPOO models in which 
the gas infall could be tuned to simulate starvation (by 
simply stopping gas infall) or ram pressure stripping (by 
removing gas already settled onto the disk). These au- 
thors concluded that ram pressure stripping is required 
to explain the truncated HI and star- forming disks of this 
galaxy. 

NGC 4736 constitutes another interesting example. At 
radii larger than 225" (~ 6 kpc), a pronounced anti- 
truncation can be clearly seen in its light profiles at 
all wavelengths. The bulge dominates the emission at 
r < 75" (^ 2 kpc) , and a prominent star-forming ring 
is visible at r ~ 40" (^ Ikpc). Both the outer disk and 
the bulge (together with the ring) were excluded when 
performing the fit, and the best- fitting model is very suc- 
cesful at reproducing the multi-band profiles of the inner 
disk, but again overestimates the true FUV and NUV 
profiles. Trujillo et al. (2009) carried out a detailed anal- 
ysis of this galaxy by fitting its sepectral energy distri- 
bution (SED) at different radii, and also by performing 
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smoothed particle hydrodynamics simulations. These 
authors favor a scenario where an oval distortion gives 
rise to both the outer disk and the enhanced central star 
formation activity, involving radial flows which of course 
the BPOO models do not consider. 

Even though the particular objects described above are 
rather extreme, Fig.lTSlshows that the model is still quite 
off in the UV for many objects with seemingly normal HI 
contents, including early-type spirals like NGC 3031 or 
NGC 2841 (see Fig. [3|). In these cases, any external in- 
teraction probably plays a minor role compared to the 
larger failings of the models for predicting gas contents. 
As mentioned at the end of Section 12. 2[ the BPOO mod- 
els parameterize the gas infall rate as a function of both 
the local mass surface density and the total mass of the 
galaxy, so that infall will proceed faster in the densest 
parts of disks and, in general, in the most massive galax- 
ies. Boissier (2000) showed that tuning the mass depen- 
dence of the infall rate can modify the present-day colors 
of galaxies. Redder stellar populations can be obtained 
if gas infall in massive early-type disks takes place even 
faster and earlier than in the finally adopted version of 
the model. Further investigation in this direction is left 
for future papers. 

6.3. Implications for the inside-out growth of disks 

After having analyzed the strenghts and weaknesses of 
the BPOO models, here we discuss the implications of the 
model fitting results regarding the inside-out growth of 
spiral disks. First, we will compare the predicted and ob- 
served color profiles in our sample, to ascertain which as- 
pects of disk evolution can be reproduced by the models 
and which ones would require a more complex approach 
(such as N-body simulations). With these limitations in 
mind, we will then study the size evolution of our disks 
since z = 1 assuming that, at least to first order, they 
have evolved as dictated by the model that best fits their 
current multi- wavelength profiles. 

6.3.1. Color profiles 

In Fig. [19] we compare the color profiles of our galax- 
ies with those predicted by the model. The gray lines 
show the color profiles within the radial range used in 
the fitting (that is, excluding the bulges). The black 
lines show the model predictions (using the KOI IMF) 
for selected values of A and Vc that roughly bracket the 
values found for most of our galaxies. The solid lines cor- 
respond to A = 0.03, and the dashed ones to A = 0.09. 
Within each pair of lines, the upper (that is, redder) 
one has Vc = 360kms~^, whereas the lower one has 
Vc — 80kms^^. Two sets of color profiles are shown: 
(FUV — 3.6 /im) profiles and {g — 3.6 /im) ones, both of 
them with and without correcting for internal extinction. 

Even after accounting for the effect of dust attenua- 
tion, the (FUV — 3.6 ^m) profiles of Sa-Sbc galaxies fall 
outside the region of the diagram delineated by the model 
predictions. This is a direct consequence of the failure of 
the models at reproducing the UV profiles in early-type 
disks, as extensively discussed in the previous sections. 
In later types, however, the agreement is excellent once 
extinction is taken into account (the difference is mini- 
mal in Sdm-Sm galaxies, since they are not particularly 
dusty) . 



The discrepancies become much less severe redward of 
the Balmer break. After correcting for internal extinc- 
tion, the model predictions for the {g — 3.6 fim) profiles 
nicely encompass the actual color profiles of our galaxies 
from early to late Hubble types. As we already showed 
in Fig. [3] with NGC 2841, the models are capable of re- 
producing the optical and near-IR profiles of early-type 
disks even when they fail in the UV. 

However, the color profiles of our galaxies become red- 
der at large galactocentric distances, thus exhibiting a 
U-like shape. Note that the fact that we can acceptably 
fit the light profiles but not the color ones is not necessar- 
ily contradictory. Our models can reproduce the globally 
exponential nature of disks, but are not sensible to small 
variations such as those that can lead to the observed 
color profiles. 

These U-shaped profiles appear to be common both in 
nearby (Bakos et al. 2008) and distant galaxies (AzzoUini 
et al. 2008a). N-body simulations show that they may 
result from a combination of a drop in the SFR (seeded 
by warps in the gaseous disk, radial distribution of an- 
gular momentum, misalignment between the rotation of 
the infalling gas and the disk, etc.) and radial stellar mi- 
gration, which would populate the outskirts of disk with 
old stars formed inwards (Roskar et al. 2008; Martinez- 
Serrano et al. 2009; Sanchez-Blazquez et al. 2009). Un- 
fortunately, none of these processes can be easily trans- 
lated into an analytic ID scheme such as ours without 
introducing too many unconstrained free variables. 

However, the simplicity of our models —compared to 
N-body simulations— comes at the advantage of being 
able to easily generate large grid of models over a wide 
range of halo masses and spins. Until similar grids of N- 
body disk simulations become available, the BPOO consti- 
tute a reasonable first-order approach to infer a galaxy's 
past evolution from its present-day photometric profiles. 

6.3.2. A first look on the past evolution of SINGS disks 

We have already checked that the model is able to 
reproduce the observed circular velocities of our galax- 
ies, as well as the expected values of the spin parameter 
—both of which seem to be unaffected by our particular 
choice of IMF. In addition, we have also carefully ex- 
plored some physical reasons that may be responsible for 
the failure of the model at reproducing the UV profiles 
of some early-type disks, which may hamper the study 
of the recent SFH in these objects. We can now proceed 
to study the evolution of the objects in our sample, by 
assuming that, at least on timescales of a few Gyr, our 
disks have evolved in a similar way as the correspond- 
ing model that best fits its present-day multiwavelength 
profiles (see also Boissier & Prantzos 2001 for a detailed 
analysis of the evolution with redshift of several physical 
properties of galaxies). 

For each model characterized with a particular pair of 
values of A and Vc, we determine the disk scale- length 
Rd at each epoch t by fitting an exponential law to the 
total stellar mass density profile: 

E,tars(r, t) = S,tars(0, t)e-'^/«^(*) (21) 

In Fig. [201 we show the temporal evolution of the disk 
scale- length for selected values of A and Vc. As expected, 
Rd increases with time in all cases. At any given epoch. 
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Fig. 19. — Comparison between the model and observed (FUV — 3.6 fira) and (g — 3.6 fira) color profiles (only those galaxies with SDSS 
data appear in the bottom panels). Gray lines show the color profiles of our galaxies, both before and after correcting for internal extinction 
("obs" and "corr", respectively). The black lines correspond to model color profiles (with the KOI IMF) for values of A and Vq that roughly 
encompass most values found in our sample. Solid lines have A = 0.03, and dashed ones have A = 0.09. Within each pair of profiles, the 
upper and redder one corresponds to Vq = 360kms~^, while the lower and bluer one has Vc = SOkms^^. 



the most extended disks are those with larger vahies of 
either parameter, as we already pointed out in Fig. [1] 

The curves describing the growth rate of i?d seem to get 
steeper with increasing Vc at fixed A. In order to quantify 
the slope of these curves, we have performed a linear fit 
to i?d as a function of time between z = 1 and z ~ Q. 
Disk growth seems to be approximately linear since z — 1 
and, as mentioned in the previous sections, it is not clear 
whether the models can describe disk evolution beyond 



that redshift, when mergers were more frequent and the 
thin disk was not fully assembled. 

In Fig. we plot the disk growth rate dRd/dt as a 
function of both the spin parameter and the rotational 
velocity. The growth rate increases up to A ~ 0.06, but 
for larger values it seems to be largely insensitive to the 
particular spin of the galaxy. It is clear that dR^/dt 
mainly depends on the circular velocity, while the spin 
parameter only seems to be relevant at large velocities. 
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Fig. 20. — Temporal evolution of the exponential scale-length of the total stellar mass density profiles. For each value of the spin 
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In general, very late-type disks appear to grow at a rate 
of 0.02-0.04kpc Gyr~^; early-type spirals, on the other 
hand, can increase their scale-lengths at a rate up to 
~ O.lkpcGyr"^, depending on their spin. 

Rather than describing the evolution of disks in terms 
of their absolute growth rate in kpcGyr"^, it is perhaps 
more illustrative to focus on their relative size increase. 
We have plotted in Fig. [22] the ratio of the scale-length of 
the stellar mass profiles at z = relative to z = 1. Inter- 
estingly, this size ratio is essentially a unique function of 
the spin, with almost no dependence on Vc —and hence 
on mass. It might seem striking that high-spin galaxies 
experience almost no change in size since z = 1. How- 
ever, as Fig. [20] demonstrates, these galaxies will already 
exhibit extended stellar mass profiles at z = 1. There- 
fore, even if the absolute growth rate is high, it will not 
have a significant impact on the relative increment in 
size. 

Can we extrapolate the conclusions obtained for our 
sample to the general population of disk-like galaxies? 
The histograms in Fig. [23] show the distribution of both 
the absolute and relative growth rates in our sample. 
Both distributions peak around the values typical for 
galaxies similar to the Milky Way, with an absolute 
growth rate of about 0.05-0.06 kpc Gyr~^ and a relative 
size increase roughly equal to 20-25% since z = 1. 

We should not blindly extend these results to the whole 
population of spiral galaxies in the Local Universe. The 
absolute growth rate depends primarily on Vc and, as 
we discussed in section [6. 1.1[ low- mass disks are consid- 
erably underrepresented in our sample. Therefore, most 



disks in a volume-limited sample would likely grow at 
slower rates than the peak value in Fig. [23] 

However, the situation is different regarding the rel- 
ative increment in size. As we stated in section 16.1.11 
the distribution of spin values in our sample matches 
reasonably well the one found in N-body simulations of 
disk formation. Therefore, we can treat our sample as 
being representative of a complete one regarding any A- 
dependent quantity. This is precisely the case of the rel- 
ative size ratio, which depends almost entirely on A alone 
according to Fig. [21] Therefore, we can safely conclude 
that most disks have probably undergone an increase of 
20-25% in their scale-lengths since z = 1 until now, re- 
gardless of their total mass. This result is in perfect 
agreement with the growth rate we estimated in Muhoz- 
Mateos et al. (2007) on a larger sample of galaxies, but 
using only extinction-corrected FUV and Ks-hand pro- 
files plus a very simple toy model of disk growth. 

It is interesting to compare our theoretical expecta- 
tions for inside-out disk growth with actual measure- 
ments of disk sizes at different redshifts. From a ob- 
servational perspective, this issue is typically addressed 
by studying the evolution (or lack of thereof) in the 
magnitude-size and mass-size relations, the average sur- 
face brightness and the size number density (Schade et 
al. 1996; Lilly et al. 1998; Simard et al. 1999; Ravin- 
dranath et al. 2004; Trujillo et al. 2004, 2006; Harden et 
al. 2005; Mcintosh et al. 2005; Trujillo & Pohlen 2005; 
Azzollini et al. 2008b). The results of such studies are 
sometimes contradictory, due to selection effects and the 
inherent difficulty in disentangling the evolution of in- 
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Fig. 21. — Growth rate of the scale-length of the stellar mass density profiles, as a function of A (a) and Vq (b). The growth rate has 
been computed by fitting as a function of t between z = \ and z = 0. Observed morphological types are coded with different symbols. 




Fig. 22. — Ratio of the stellar disk scale-lengths at 2: = and 2: = 1 as a function of the spin (a) and the rotational velocity (b). Observed 
morphological types are coded with different symbols. 



dividual galaxies from the evolution of a population of 
galaxies as a whole. 

Boissier & Prantzos (2001) confronted their predicted 
size-luminosity trend with several observed data-sets 
from the literature. Here we revisit this issue, and com- 
pare our model with the mass-size relation derived by 
Barden et al. (2005). These authors determined disk ef- 
fective radii at various z by fitting a Sersic model (Sersic 
1968) to HST images from GEMS (Galaxy Evolution 
from Morphologies and SEDs; Rix et al. 2004). Stel- 
lar masses were derived from SED fitting to COMBO-17 



data for the same objects. As a local anchor for their 
study they relied on a sample of nearby galaxies from 
SDSS. These authors found little or no evolution with 
z in the mass-size relation. Assuming that galaxies can 
only become more massive with time, they argued that 
they should also increase their sizes accordingly. 

To replicate these measurements, for each one of our 
model disks we fitted a Sersic profile to the radial distri- 
bution of stellar mass at different epochs. The results are 
shown in Fig. [24l Each track follows the evolution of the 
effective radius and total stellar mass of a model galaxy 
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Fig. 24. — Evolution of the effective radius and total stellar mass 
of our model galaxies. The effective radius was obtained by fitting 
a Sersic model to the stellar mass profiles at different z. Each 
track follows the evolution of a model disk from z = 1 to z = 
(the 2 = step is marked with a symbol). The corresponding spins 
and circular velocities of each model are indicated with different 
line styles and symbols, respectively. The closed curve shows the 
observed mass-size relation from Barden et al. (2005). 



from 2 = 1 to z = (the z = step is marked with a 
symbol). Models with different velocities and spins are 
coded with different symbols and line styles, respectively. 
The irregular closed line encompasses the empirical data- 
points of Barden et al. (2005), including different redshift 
bins between z = and z = 1 (see their Fig. 10). 

Our model is in perfect agreement with the observed 
mass-size trend. Note that the effective radii in Barden 



et al. (2005) were corrected to the restframe V band, 
without any further correction to get the actual stellar 
size. Sersic fits to our model profiles indicate that the 
V-band effective radii are typically 0.05-0.15 dex larger 
than the stellar ones, so the closed curve in Fig.lMlshould 
be shifted downwards by that amount. The agreement 
would be still excellent; it would actually improve if we 
note that most galaxies are expected to lie between the 
A = 0.02 and A — 0.05 tracks, according to the A proba- 
bility distribution shown in Fig. 4. 

7. CONCLUSIONS 

In this paper we have fitted the UV, optical and near- 
IR profiles of 42 disk-like galaxies of the SINGS sample 
with the models of Boissier & Prantzos (1999, 2000). 
In order to recover the unattenuated starlight profiles 
at each wavelength, we have used the radially-varying 
TIR/UV ratio as tracer of the internal extinction. The 
disk evolution models are calibrated on the Milky Way 
(BP99) and further extended to other disk-like galaxies 
through scaling laws derived from the ACDM scenario 
(BPOO). By considering in a consistent framework the 
gas infall, the star formation activity and the subsequent 
chemical evolution, the models are capable of predicting 
the current multi-band profiles of spirals as a function of 
only two parameters: the maximum circular velocity of 
the rotation curve, Vc, and the dimensionless spin pa- 
rameter, A. By fitting the profiles of the SINGS galaxies 
with those predicted by the models, we have not only 
checked the accuracy of the models themselves, but have 
also used them to infer the size evolution of the SINGS 
galaxies. The main conclusions of this work can be sum- 
marized as follows: 

1. The rotational velocities are in good agreement 
with those measured from observed HI rotation 
curves, as well as with those estimated from the 
TuUy-Fisher relation. In the latter case, the val- 
ues of Vc derived from the model fitting tend to 
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be 10-20% larger than those predicted by the TF 
relation, but mainly for the most massive disks. 

Most galaxies in our sample exhibit spin values of 
A ^ 0.03. In fact, even though the sample is not 
volume-limited, its statistical distribution of spin 
values closely resembles the narrow distributions 
usually found in N-body simulations, which typi- 
cally have an almost universal peak at A ^ 0.03- 
0.04. 

There is a clear, well-known trend between Vq and 
Hubble type, in the sense that early-type disks 
have larger circular velocities —and are hence more 
massive— than late-type ones. There is not, how- 
ever, any evident trend between the morphological 
type and A, which supports the findings of numer- 
ical simulations that most haloes posses the same 
spin, regardless of their total mass or mass assem- 
bly history. 

While there is an excellent agreement between the 
model predictions and the observed profiles in the 
optical and near-IR bands, significant departures 
may arise in the UV bands, depending on the mor- 
phological type and the particular choice of IMF. 
The Kroupa (2001) IMF yields excellent results in 
Sc-Sd spirals, but overestimates the UV luminos- 
ity in early- type disks, and to a much lesser extent 
in Sdm-Sm ones. The Kroupa et al. (1993) IMF 
brings the UV model profiles into better agreement 
with the observed ones in Sb-Sbc spirals, as well as 
in Sdm-Sm ones, but at the expense of loosing the 
excellent fits for Sc-Sd disks. While differences in 
the high-mass end of the IMF might indeed play a 
role in very late-type galaxies, it is doubtful that 
the IMF is behind the discrepancies in the UV pre- 
dictions for early-type disks. Anyway, the values of 
A and Vc are largely unaffected by the specific IMF 
chosen. For the Hl-deficient galaxies in our sample, 
gas removal due to interactions with the intraclus- 
ter medium is the most likely culprit. For those 
galaxies with normal HI masses it may be neces- 
sary to revisit the mass dependence of the gas infall 
rate, since the model seems to retain too much gas 
in these objects. 

The metallicity gradients predicted by the mod- 
els are ~ 0.015 dexkpc"^ steeper than the observed 
ones. The central oxygen abundances depend on 
the IMF: the values yielded by the K93 IMF are in 
perfect agreement with the observed central metal- 
licities, but those obtained with the KOI one over- 
estimate the real values by ~ 0.7 dex. 

According to the models, the absolute growth rate 
(in kpcGyr"^) of the exponential scale-length of 
disks depends mainly on Vc, with rapidly rotating 
disks expanding faster. In our sample, most galax- 
ies have their scale-lengths increased by about 0.05- 
0.06 kpc each Gyr. Still, this is not representative 
of the overall population of disks, since low-mass 
ones are underrepresented in our sample. 

The ratio between the current disk scale-length and 
that at z = 1 is a decreasing function of A, with 



little dependence on Vq. Even though high-spin 
disks grow faster in absolute terms, such a rapid ra- 
dial expansion does not significantly increase their 
scale-lengths, which are already considerably large 
at z = 1. On average, most disks in our sample 
are now 20-25% larger than at z = 1. This value 
can be treated as being representative of a volume- 
limited sample, given that our galaxies have the A 
distribution expected for such a kind of sample. 

8. The model predicts that disk galaxies should simul- 
taneously increase their sizes and stellar masses as 
time goes by. The results of the model for a grid 
of values of A and Vc provide a perfect match to 
the observed constancy of the mass-size relation 
between z — and z — 1. 
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APPENDIX 

TWO-DIMENSIONAL DISTRIBUTION OF VALUES 

When fitting multi-wavelength profiles with the disk evolution models of BPOO, one must bear in mind that the 
circular velocity and the spin may not act as completely independent parameters. Depending on the particular shape 
of a galaxy's profile at different wavelengths, variations in one parameter might be compensated by variations in the 
other one while still providing an acceptably good fit. 

In order to depict the possible internal degeneracies between A and Vq, in this appendix we present the two- 
dimensional distributions obtained for each galaxy in our sample (Fig. I25|) . Even though we keep track of the 
individual distributions corresponding to each particular band for each galaxy, the plots presented here show the 
distribution of total values taking into account all bands. The best fitting model is marked with a white dot, and 
the area encompassing all models with x ^ ^Xmin has been delimited with a white closed line. The same range in A 
and Vq is displayed in all panels, except for the small subset of galaxies requiring larger spins and/or lower circular 
velocities than those in our initial grid. For the sake of clarity, those galaxies have been grouped together at the end. 

In general, some galaxies exhibit some degree of anticorrelation between both parameters for low and intermediate 
values of A, since the increment in the radial scale- length caused by augmenting A can be partly compensated by 
decreasing Vc- In some other objects, the x^ distribution around the best fitting model does not show any significant 
degeneracy. Finally, for large values of the spin the correlation is possitive: further incrementing A significantly 
decreases the central surface brightness, which can be compensated to some extent by increasing Vc —even though 
this tends to augment the radial scale-length as well. 

The fact that the spin parameter is not as strongly constrained as the circular velocity is mainly due to the different 
way in which both quantities affect the radial profiles (Fig. 1). Modifying Vc will shift the model profiles above or 
below the observed ones, thus rapidly increasing the x^ value. Varying A, on the other hand, will mainly change the 
scale-length alone. Given that the observed profiles exhibit inhomogeneities that deviate from our smooth predictions, 
this leaves some room for varying A while still obtaining a good fit. 
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TABLE 2 
Model results. 
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Note. — Results from the model fitting. (1): Galaxy name. (2), (6): Dimensionless spin parameter. (3), (7): Maximum circular velocity. 
(4), (8): Central value of 12 + \og{0/H). (5), (9): Radial metallicity gradient. (10): Temporal growth rate of the stellar disk scale-length, 
obtained by fitting Rd{t) between z = 1 and z = 0. (11): Ratio of the stellar disk scale-lengths at z = and z = 1. Neither (10) nor (11) vary 
noticeably with the IMF. 
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